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1 Introduction 

Understanding the electron dynamics and transport in metallic and semicon- 
ductor nanostructures - such as metallic nanoparticles, thin films, quantum 
wells and quantum dots - represents a considerable challenge for today's con- 
densed matter physics, both fundamental and applied. 

Experimentally, thanks to the recent development of ultrafast spectroscopy 
techniques, it is now possible to monitor the femtosecond dynamics of an 
electron gas confined in metallic nanostructures such as thin films [1, 2, 3, 4, 
5, 6, 7, 8], nanotubes [9], metal clusters [10, 11] and nanoparticles [6, 7, 12, 13]. 
Therefore, meaningful comparisons between experimental measurements and 
numerical simulations based on microscopic theories are becoming possible. 

The dynamics of an electron gas confined in a metallic nanostructure is 
characterised by the presence of collective oscillations (surface plasmon) whose 
spectral properties depend on several conditions of temperature, density, and 
coupling to the environment. At lowest order, the linear response of the elec- 
tron gas is simply given by the plasma frequency ujp = (e^n/meo)^/^ (up to 
a dimensionless geometrical factor), and does not depend on the temperature 
or the size of the nano-object. The plasma frequency represents the typical 
oscillation frequency for electrons immersed in a neutralizing background of 
positive ions, which is supposed to be motionless because of the large ion mass. 
The oscillations arise from the fact that, when some electrons are displaced 
(thus creating a net positive charge), the resulting Coulomb force tends to 
pull back the electrons towards the excess positive charge. Due to their in- 
ertia, the electrons will not simply replenish the positive region, but travel 
further away thus re-creating an excess positive charge. This effect gives rise 
to coherent oscillations at the plasma frequency. Notice that, for a metallic 
nanostructure, the inverse plasma frequency is typically of the order of the 
femtosecond - this coherent regime can therefore be explored with the ultra- 
fast spectroscopy techniques developed in the last two decades. 
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The coherence of such collective motions is progressively destroyed by Lan- 
dau damping (i.e. by coupling to the internal degrees of freedom of the electron 
gas) and by electron-electron or electron-phonon collisions. The damping of 
the plasmon was observed experimentally in gold nanoparticles [14] and was 
studied theoretically in several works [15, 16, 17]. 

Although the linear response of the surface plasmon has been known for a 
long time, fully nonlinear studies have only been performed in the last decade 
and have revealed some interesting features. Our own contribution to this 
research area has mainly focussed on the nonlinear electron dynamics in thin 
metal films, where the emergence of ballistic low-frequency oscillations has 
been pointed out [18]. 

On the other hand, the same type of collective electron motion is also 
observed in semiconductor nanostructures, such as quantum wells and quan- 
tum dots. Although the spatial and temporal scales differ by several order of 
magnitudes with respect to metallic nanostructures (due the large difference 
in the electron density), the relevant dimensionless parameters take similar 
values in both cases [19]. For instance, the effective Wigncr-Seitz radius is of 
order unity for both metallic and semiconductor nano-objects. Therefore, the 
electron dynamics can be investigated using similar models and both types of 
nano-objects are expected to share a number of similar dynamical properties. 

In this review article, we will describe the collective electron dynamics in 
metallic and semiconductor nanostructures using different, but complemen- 
tary, approaches. For small excitations (linear regime), the spectral properties 
can be investigated via quantum mean-field models of the TDLDA type (time- 
dependent local density approximation), generalized to account for a finite 
electron temperature. In order to explore the nonlinear regime (strong exci- 
tations), we will adopt a phase-space approach that relies on the resolution of 
kinetic equations in the classical phase space (Vlasov and Wigner equations). 
The phase-space approach provides a useful link between the classical and 
quantum dynamics and is well suited to model effects beyond the mean field 
approximation (electron-electron and electron-phonon collisions). We will also 
develop a quantum hydrodynamic model, based on velocity moments of the 
corresponding Wigner distribution function: this approach should lead to con- 
siderable gains in computing time in comparison with simulations based on 
conventional methods, such as density fimctional theory (DFT). 

The above studies all refer to the charge dynamics in a semiconductor or 
metallic nanostructure, which has been intensively studied in the last three 
decades. In more recent years, there has been a siirgc of interest in the spin 
dynamics of the carriers, mainly for possible applications to the emerging 
field of quantum computing [20]. A promising approach to the development 
of a quantum computer relies on small semiconductor devices, such as quan- 
tum dots and quantum wells [21]. To implement basic qubit operations, most 
proposed schemes make use of the electron spin states, so that a thorough 
understanding of the spin dynamics is a necessary prerequisite. Nevertheless, 
in order to manipulate the electrons themselves, one must necessarily resort to 
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electromagnetic fields, which in turn excite the Coulomb mean field [22, 23]. 

The charge and spin dynamics arc therefore closely intertwined and both must 
be taken into account for a realistic modelling of semiconductor-based qubit 
operations. 

The ultrafast magnetization (spin) dynamics in ferromagnetic nanostruc- 
tures has also attracted considerable experimental attention in the last decade. 
Pioneering experiments [24] on ferromagnetic thin films revealed that the mag- 
netization experiences a rapid drop (on a femtosecond time scale) when the 
films is irradiated with an ultrafast laser pulse, after which it slowly regains 
its original value on a time scale close to that of the electron-phonon cou- 
pling. Despite many attempts [24, 25, 26], a clear theoretical explanation for 
these effects is still lacking. Here, we will illustrate how this problem can be 
addressed using some of the techniques developed for the electron dynam- 
ics, particularly quantum mean-field and phase-space methods, which will be 
generalized to include the spin degrees of freedom. 

2 Models for the electron dynamics 

Metallic and semiconductor nano-objects operate in very different regimes, 
as the electron density is several orders of magnitudes larger for the former. 
Consequently, the typical time, space, and energy scales can be very different, 
as illustrated in Table 1. However, if one takes into account the effective 
electron mass and dielectric constant, the relevant dimensionless parameters 
turn out to be rather similar [19]: for instance, from Table 1 we see immediately 
that the ratio of the screening length (iscreen = vf l^p, where vf is the Fermi 
velocity) to the effective Bohr radius as = ^Treh^/me'^ is of order unity. The 
same happens for the ratio of the plasmon energy hujp to the Fermi energy 
Ef, so that the normalized Wigner-Seitz radius is also of order unity for 
both cases. ^ 

It is therefore not surprising that the electron dynamics of both types of 
nanostructures can be described by means of similar models. A bird's-eye view 
of the various relevant models is provided in Fig. 1. The diagram represents 
the various levels of modeling for the electron dynamics, both quantum (left 
column, orange) and classical (right column, blue). The highest level of de- 
scription is the A''-body model, which involves the resolution of the A''-particle 
Schrodinger ccjuation in the quantum regime, or the TV-particle Liouville equa- 
tion for classical problems (the latter is of course equivalent to Newton's equa- 
tions of motion). This is a difficult task even classically, although molecular dy- 
namics simulations that solve the exact A''-body problem can nowadays attain 



^ For a quantum well, all relevant lengths far exceed the semiconductor lattice spac- 
ing aiatt — 5 A. This makes semiconductor systems a much better approximation 
to jellium (i.e., a continuum ionic density profile) than simple metals, for which 
the lattice spacing is comparable to the other electronic lengths. 
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Table 1. Typical time, space, and energy scales for metallic and semiconductor 
nanostructures 
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a considerable level of sophistication. For Newton's equations with two-body 
interactions, the numerical complexity grows at most as iV^, and in some cases 
this can be reduced to a logarithmic dependence. Quantum-mechanically, the 
A''-body problem is virtually unmanageable, except for very small systems, 
because the size of the relevant Hilbert space grows exponentially with N . 
Nevertheless, exact simulations of the A''-body Schrodinger equation can be 
performed using the so-called configuration interaction (CI) method. We have 
used this approach to study the exact electron dynamics in semiconductor 
quantum dots containing up to four electrons. 

For larger systems, some rather drastic approximations need to be made 
if we want to end up with a mathematically and numerically tractable model. 
Most such reduced models are improvements on the so-called 'mean field ap- 
proximation', which states that the motion of a single electron is determined 
by the positions and velocities of all other particles in the system. Such collec- 
tive behavior is possible because of the long-range nature of electromagnetic 
forces. The mean field approach can be viewed as a zeroth-order approxima- 
tion to the iV-body problem in which two-body (and higher order) correla- 
tions between the particles have been neglected. Classically, this procedure is 
known as the BBGKY hierarchy (from the names of Bogoliubov, Born, Green, 
Kirkwood, and Yvon) [27]. 

For classical systems of charged particles (plasmas), the mean-field dy- 
namics is governed by the Vlasov equation, which describes the evolution of a 
one-particle probability density in the phase space. The quantum analog of the 
Vlasov equation is provided by the time-dependent Hartree equations, which 
are actually one-body Schrodinger equations evolving in the mean-field poten- 
tial. In both cases, the mean field is obtained by solving Maxwell's equations, 
often reduced, in the electrostatic limit, to the sole Poisson's equation. 

In this review, we concentrate on quantum mechanical models. Several im- 
provements have been proposed to the Hartree equations (which were derived 
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in 1927, just one year after Schrodinger's seminal paper on the wave equation), 
most notably Feck's correction (1930). Indeed, the Hartrcc method does not 
respect the principle of antisymmetry of the wavefunction, although it does 
use the Pauli exclusion principle in its less stringent formulation, forbidding 
the presence of two electrons in the same quantum state. The Hartrcc-Fock 
equations respect the antisymmetry of the wavefunctions, thus leading to an 
extra interaction term between the electrons, termed the 'exchange interac- 
tion'. 

A particularly successful extension of the mean-field approach is the 
density- functional theory (DFT), which was developed by Hohenberg, Kohn, 
and Sham in the mid 1960s [28]. Originally developed for the gromid state 
at zero temperature, it has subsequently been extended to finite temperature 
and time-dependent problems. As its name suggests, DFT states that all the 
properties of a many-electron systems are determined by the electron spatial 
density, rather than by the wavefunctions. DFT allows one to introduce in 
the mean-field formalism effects that go beyond the strict mean-field approx- 
imation, particularly the exchange interaction described above. Indeed. DFT 
can deal with higher order correlations between the electrons, in principle ex- 
actly if the exact density functional were known. In practice, one has to make 
an educated guess for the appropriate correlation functional, which leads to 
various empirical approximations. Nevertheless, DFT has proven immensely 
useful for a wide range of electronic structure calculations. 

The Hartrce equations can be cquivalcntly recast in a phase-space formal- 
ism by making use of the Wigner transformation, which was introduced by 
E. Wigner in 1932 [29]. The resulting Wigner function is a pseudo probability 
distribution, which can be used to compute expectation values just like its 
classical counterpart. Unfortunately, the Wigner function can take negative 
values, which precludes the possibility of interpreting it as a true probability 
density. 

By taking velocity moments of the Wigner equation - and using some ap- 
propriate closure hypotheses - one can derive a set of quantum hydrodynam- 
ical (or fluid) equations that govern the evolution of macroscopic quantities 
such as the particle density, average velocity, pressure, heat flux etc. Com- 
pared to the Wigner approach, the hydrodynamical one is obviously numeri- 
cally advantageous, as it requires the resolution of a small number of equations 
in real (not phase) space. Generally speaking, hydrodynamical methods yield 
accurate results over distances that are larger than the typical electrostatic 
screening length, which is the Debye length Ad = {ksTee / e^nY/"^ for classical 
plasmas and the Thomas- Fermi screening length Lp = vp/i^p for degenerate 
electron gases (see Table 1). 

In the following subsections, we shall present a brief overview of most of 
the quantum models featuring in Fig. 1. 
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Fig. 1. Bird's-eye view of the models used to describe the electron dynamics. From 
top to bottom: A^'-body, mean-field, and macroscopic (hydrodynamic) theory. Left 
column (orange): quantum models; right column (blue): classical models. Notation: x 
= exchange; xc = exchange and correlations; \d = Debye length (classical screening 
length); Lf = Thomas- Fermi screening length; k = typical wavevector; BBGKY = 
Bogoliubov, Born, Green, Kirkwood, Yvon hierarchy. 
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2.1 Exact iV-body simulations: the Configuration Interaction (CI) 
method 

Method 

In the Hartree-Fock model (HF), the many-body wave- function is approxi- 
mated by a single Slater determinant leading to a correlation between electrons 
having the same spin. However, electrons of different spin are not correlated in 
this approximation. This is why the difference between the exact value of the 
energy, and the HF value is called the correlation energy. There are a number 
of quantum chemistry methods, which attempt to improve the description of 
the many-body wave-function. The most important one is the so-called con- 
figuration interaction method (CI) [30] which is based on the minimization 
of the energy with respect to the expansion coefficients of a trial many-body 
wave- function expressed as a linear combination of Slater determinants. With 
respect to the models based on density functional methods the drawback of 
the CI method is its unfavorable scaling with the system size. Indeed, the di- 
mension of a full CI calculation grows factorially with the number of electrons 
and basis functions. 

From the above considerations, it is clear that CI calculations are restricted 
to confined systems with very few electrons (typically less than 10). In quan- 
tum chemistry, the "basis set" usually refers to the set of (nonorthogonal) 
one-particle functions used to build molecular orbitals. Concerning the com- 
putational methodology for confined electron systems (atoms, molecules, clus- 
ters, nanoparticles, quantum dots...) localized basis sets arc the traditional 
choice and the most common type of basis functions is the Gaussian func- 
tions. It is worth noticing that, from the knowledge of the exact many-body 
wave-function, one can in principle: (i) compute the temporal evolution of 
the system, including the dynamical correlations; (ii) obtain the true excited 
states of the system. 

In the following, an application of the CI method in the field of semicon- 
ductor nanostructures and quantum dots is presented. 

Application 

Recent progress in semiconductors technology allows the realization of quan- 
tum systems composed of a small number of electrons (even a single electron!) 
confined in nanometer-scale potential wells. These systems, which provide 
highly tunable structures for trapping and manipulating individual electrons, 
are often named artificial atoms or quantum dots and arc good candidates 
for the emerging technology of quantum computing. They have certain simi- 
larities with atoms in the sense that they have a discrete electronic structme 
that follows the well-known Hund's rule of atomic physics. However, in quan- 
tum dots the electrons are generally confined by harmonic or quasi-harmonic 
potentials, whereas atoms are characterized by Coulomb confinement poten- 
tials. The spectral properties of quantum dots are exotic with respect to the 
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properties of atoms in the sense that most of the oscillator strength is concen- 
trated almost exclusively on one dipolar transition. This property is a direct 
consequence of Kohn's theorem (KT) and does not depend on the number of 
electrons, the strength of the confinement or the electron-electron interaction 
[31]. 

In a recent work [32] , we have investigated quasi- two-dimensional Gaussian 
quantum dots containing up to four electrons within the framework of the CI 
method which allows in principle an exact treatment of the many-electron 
system. The Schrodinger equation for A''-electrons confined by a potential 
Vext is given by 

H^il,...,N) = E^il,...,N) (1) 

where (1, iV) represents the space [r, = {xi,yi, Zi)] and spin coordinates of 
the electrons and 

« = i:-|sV.' + E3i;^ + EWr.). (2) 

1=1 i>j ' ■' ' 2=1 

The confinement is modelled by an external one-particle anisotropic Gaussian 
potential given by 

V^xtivi) = -£>exp [-7(ar? + yf)] + ^mW^zf . (3) 

It is worth noticing that for suflaciently large values of Uz the electrons of 
the system are strongly compressed along the z direction. Therefore, in this 
situation, the system can be regarded as a quantum system confined by a two- 
dimensional Gaussian- type potential, i.e., as a quasi-two-dimensional Gaus- 
sian quantum dot. Since a Gaussian potential can be approximated close to 
its minimum by an harmonic potential, the potential of Eq. (3) is suitable 
for the modelling of anharmonic quantum dots. The anharmonicity of the 
confinement can be characterized by the depth of the Gaussian potential D 
and by the quantity u) = ^J2D^/m. Thus, when D is much larger than hu) 
the Gaussian potential has many bound states and the potential curve follows 
closely the harmonic oscillator potential leading to a small anharmonicity of 
the system. On the other hand, when D is slightly larger than huj the Gaus- 
sian potential has only few bound states and, therefore, deviates strongly from 
the harmonic potential leading to a large anharmonicity. Also, a large (small) 
value of io corresponds to a strong (weak) confinement with respect to the 
electron-electron interaction. 

The wave-function is approximated by a linear combination of cartesian 
anisotropic Gaussian- Type Orbitals (c-aniGTO) [33]. A c-aniGTO centered 
at {bx, by,bz) is defined as 

X"'^(r,6) = xiy,^^z^,: exp(-C.< - - (4) 

where Xb^ = {x — bx) etc... Following the quantum chemical convention the 

orbitals are classified as s-type, p-type, for I = ax + (ly + = 0, 1, respec- 
tively (this sum controls the value of the orbital angular momentum). The 
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{bx,by,bz) parameters have been chosen to coincide with the center of the 
confining potential. This type of basis sets was found to be the most suitable 
one for expanding the eigenfunctions of an electron in an anisotropic harmonic 
oscillator potential. The calculations have been performed using the OpenMol 

Program ^. 

Energy spectra and oscillator strengths have been calculated for different 
strength of confinement oj and potential depth D. The effect of the electron- 
electron interaction on the distribution of oscillator strengths and the break- 
down of the KT has been examined by focusing on the results with the same 
value of D/hui i.e. with the same anharmocity. 

A substantial red-shift has been observed for the dipolc transitions cor- 
responding to the excitation into the center-of-mass mode. The oscillator 
strengths, which are concentrated exclusively in the center-of-mass excita- 
tion in the harmonic limit, arc distributed among the near- lying transitions 
as a result of the breakdown of the Kohn's theorem. The distribution of the 
oscillator strengths is limited to the transitions located in the lower-energy 
region when lu is large (i.e. for strongly confined electrons) but it extends to- 
wards the higher-energy region when u becomes small (i.e. for weakly confined 
electrons). 

The analysis of the CI wave functions shows that all states can be classified 
according to the polyad quantum number Vp [32]. The distribution of the 
oscillator strengths for large ui occurs among transitions involving excited 
states with the same value of Vp as the centcr-of-mass excited state, Vp^cmi 
while it occurs among transitions involving the excited states with Vp — '^p^cm 
and Vp = Vp^cm+2 for small u). 

2.2 Time-Dependent Density Functional Theory (TDDFT) and 
the Local-Density Approximation (LDA) 

Time dependent density functional theory (TDDFT) extends the basic ideas 
of static density functional theory (DFT) to the more general situation of 
systems under the influence of time dependent external fields. Tliis dynamical 
approach relies on the electron density n{r,t) rather than on the many-body 
wave function !Z'(ri, r2, rjv, i) of the system. In fact, the central theorem 
of the; TDDFT is the Runge-Gross theorem [34, 35, 36] which tells us that all 
observables are uniquely determined by the density. 

Prom the computational point of view, with respect to the resolution of the 
time dependent Sc;lirodinger equation (TDSE) of an A^-electron system, the 
complexity is strongly reduced when using TDDFT. Indeed, the wave func- 
tion depends on 3A'' -|- 1 variables (ri, r2, rjv, t) while the density depends 
only on 4 variables {r,t). This is one of the reasons why this method has 
become so popular. A practical scheme for computing n{r,t) is provided by 



* see http://www.csc.fi/gopenmol 
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the Kohn-Sham (KS) formulation of the TDDFT [28]. In the latter, noninter- 

acting electrons are moving in an effective local potential constructed in such 
a way that the KS density is the same as the one of the interacting electron 
system. The advantage of this formulation lies in its computational simplic- 
ity compared to other quantum-chemical methods such as time-dependent 
Hartree-Fock or configuration interaction. The KS equations read as 

in^^<j)k{r,t) = (^-^V^ + Veff{r,t)^ <f)k{r,t) (5) 

with the KS density 

n{r,t)=Y,fk\Mr,t)\' (6) 
fe=i 

where fk denotes the occupation numbers of the ground state, and 

Veff{r,t) = Ve^t{r,t) + VH{r,t) + VUr,t) . (7) 

In the above expression the first term is the external potential (ionic potential, 
laser field...), the second is the Hartree potential, which is a solution of the 

Poisson's equation, and the last term is the exchange-correlation potential. 

The most popular choice for V^c is the so-called adiabatic local density 
approximation (ALDA) given by 

yxc{r,t) = ^ [nexc(n)]„=„(r,t) , (8) 

where exc{n) is the exchange-correlation energy density for an homogeneous 
electron gas of density n. In this approach, the same functional used to cal- 
culate the properties of the ground state is employed in the dynamical simu- 
lations. 

The validity of the local approximation has been discussed in many papers 
and textbooks [37] . This approximation works remarkably well for inhomoge- 
neous electron systems. In contrast, the validity of the adiabatic approxima- 
tion has been less thoroughly analyzed. Generally speaking, this approach is 
expected to hold for finite systems and for processes that evolve very slowly 
in time. The situation in bulk solids is more controversial since significant 
deficiencies in the description of absorption spectra have been noticed [38] . It 
was shown by Dobson [39] that ALDA fulfills the Kohn theorem when applied 
to a system of interacting electrons confined in an external parabolic poten- 
tial. This theorem guarantees the existence of a collective state at the same 
frequency as the harmonic potential. It corresponds to a rigid oscillation of 
the many-body wavefunction around the center of the external potential. 

Only a few attempts have been made to go beyond ALDA. To date, the 
most important ones are the work of Gross and Kohn [40] and that of Vignale 
and Kohn [41], the latter being the most promising in particular for studying 
electron relaxation phenomena [42]. Contrarily to ALDA, the approach of 
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Gross and Kohn, which uses a frequency-dependent parametrization of the 
exchange-correlation kernel (sec below), does not fulfill the Kohn theorem 
[31, 39]. This problem was further investigated by Vignale and Kohn [41], 
who proposed a new theory based on the so-called current density functional 
theory (CDFT). This model is described in detail in [43]. CDFT was originally 
derived by Vignale and Rasolt [44] to describe, within the framework of DFT, 
situations where strong magnetic fields and orbital currents cannot be ignored. 

Few works have been devoted to the study of the nonlinear electron dy- 
namics in finite metallic systems exposed to strong external fields. Indeed, the 
resolution of the time-dependent Kohn-Sham equations (5) is a very difficult 
task particularly for 3D systems. Some pioneering work on free simple-metal 
clusters was performed by E. Suraud in Toulouse and P.-G. Reinhard in Erlan- 
gen [45]. More recently, Gervais et al. [46] have investigated the same problem 
in 3D geometry using a spherical basis expansion technique. This approach is 
restricted to small metal clusters. The interaction of strong femtosecond laser 
pulses with a Ceo molecule (which possesses 240 delocalized electrons and can 
therefore be considered as a metallic nano-objcct [47]) has been investigated 
in Ref. [48] by employing a TDDFT approach. Still concerning the fuUerene 
molecule, Cormier et al. [49] studied multiphoton absorption processes by solv- 
ing numerically the associated time-dependent Schrodinger equation (TDSE) 
in the single active electron (SAE) approximation. This approximation con- 
sists in solving the equations (5) by using, instead of the time-dependent 
effective potential Veff{r,t) given in Eq. (7), the static effective potential of 
the ground state together with the time-dependent electric potential of the 
laser. 

Let us now examine the linear regime, which has received much wider 
attention in the past. 

Under the condition that the external field is weak, the simplest way to 
implement TDDFT is to work within the framework of the linear response 
theory. This approximation was first introduced by Zangwill and Soven [50] 
in the context of atomic physics for the study of photoionization in rare gases. 
Subsequently, this formalism has been successfully extended to the study of 
more and more complex electron systems: molecules [51], simple metal clusters 
[52], noble metal clusters [53], thin metal films [54], quantum dots [55], and 
condensed phase systems [38]. 

To date, in the field of nanoparticle physics, most applications of the time- 
dependent Kohn-Sham formalism have been performed at zero electron tem- 
perature. In order to interpret time-resolved pump-probe experiments carried 
out on noble metal nanoparticles, we have recently extended this approach to 
finite temperature. In the following we provide a brief overview of the model 
with the basic equations. 
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Ground state 

The electron gas is assumed to be at thermal equilibrium with temperature 
Te. In the Kohn-Sham formulation of the density functional theory at finite 
temperature within the grand-canonical ensemble [56] , the ground-state elec- 
tron density n of an A^-electron system is written, in terms of single-particle 
orbitals and energies Si, as 

oo oo 

n{r) = ^ /fc nfc(r) =J2h \Mr)\^ (9) 
fe=i fe=i 

where fk = [1 + exp{{ek — lJ.)/kBTe}]~^ are the Fermi occupation num- 
bers and n is the chemical potential. These orbitals and energies obey the 
Schrodinger equation 



Mr)=eiMr), (10) 



where 14// (r) is an effective single-particle potential given by 

Veffir) = Ve^tir) + Vnir) + V^,{r) , (11) 

where Vextif) is an external potential (e.g. due to the ionic background), 
Vff(r) is the Hartree potential solution of Poisson's equation, and Vxc{r) is 
the exchange-correlation potential defined by 

^"""^^^ " ["'^^c(")]„=„(r) , (12) 

where Qxcin) = / n{r) LUxc{n{r)) dr is the exchange-correlation thermody- 
namic potential [57] . The temperature appears in the self-consistent procedure 
only through the occupation numbers and the exchange-correlation thermo- 
dynamic potential. 

For low temperature (i.e. Tg < TF[n{r)] where rF[n(r)] = (37r^n(r))^^^ 
is the local Fermi temperature), u)xc{n) may be safely replaced by its value at 
Te = 0, i.e. by exc{n). The chemical potential is determined sclf-consistently 
by requiring the conservation of the total number of electrons from Eq. (9) 
[58, 59]. 



Excited States 

In the usual first-order TDLDA at Te = in the frequency domain, the in- 
duced electron density Sn{r; u) is related to SVextir'; lo), the Fourier transform 
(with respect to time) of the external time-dependent potential (generated, 
for instance, by the electric field of a laser beam), via the relation [50, 52, 60] 
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where xi'''^ '^'j ^) is the retarded density correlation function or the dynamic 
response function. It is possible to rewrite the induced density as 

5n{r;ij) = j x°{r,r';io) SVeff{r';ij) dr' (14) 

with 

5V^ff{r;.) = SVMr;^:) + ^-^^ dr' 

+ J f.c{r,r';uj) 6n{r';uj) dr' (15) 

where the function fxc{f", r' \ uj) is the Fourier transform of the timc-dcpcndent 
kernel defined by fxc{j',t;r',t') = 5Vxc{r,t)/6n{r',t') and x°('",'"';w) is the 
non-interacting retarded density correlation function. Prom Eqs. (13)-(15) we 
see that x° and x are related by an integral equation (Dyson-type equation) 



X(r, r'; uj) = x°(r, r'; ^) + J J X {r, r"; oj) 

X K{r",r"';uj) x{r"',r';uj) dr"dr"', (16) 
with the residual interaction defined by 

i^(r,r';u.) = ^— ^ + /.e(r,r';u;). (17) 

In the adiabafAc local-density approximation (ALDA) the exchange-correlation 
kernel is frequency-independent and local, and reduces to [50, 60] 

Uc{r, r') = ^ [Kc(n)]_„(,) S (r - r') . (18) 

It should be mentioned that the functional, V^c in the above equation is the 
same as the one used in the calculation of the ground state [see Eq. (12)]. For 
spin-saturated electronic systems, we have 



X 



V,r';a;) = 2^ [/° - /°] 



<P*{r)Mr)rk{r')<l>j{r') 



JK 



J2 ^Ur)Mr') G+ir,r';ek + tKo) + 

k 

occ 

J2Mr)(t>Ur')Gl{r,r';ek-nuj) (19) 
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where (j)k (r) and £fc are the one-electron Kohn-Sham wave functions and ener- 
gies, respectively. G_|_ is the onc-particlc retarded Green's function and are 
the Fermi occupation numbers at = K (0 or 1). All the above quantities 
are obtained with the procedure described in the preceding subsection with 
fk — fk Eq. (9). In order to produce numerically tractable results, wc have 
added a small imaginary part to the probe frequency, so that ui ^ ui + iS with 
T] = %5. 

At finite electron temperature, the grand-canonical non-interacting re- 
tarded density correlation function reads [61] 

X°(r, r'; a;; Te) = ^ V exp | - -1^ [E^{N) - Nf,] ] 

xx%{r,r'-u;;T„) (20) 
where Zq is the grand-canonical partition function 

= ^ exp I - ^ [E„ {N) - N^^] | (21) 

with En{N) the energy of the state \nN) having A'' electrons, the chemical 
potential and 

/ / {nN \h{r)\ mN) {mN \hir')\nN) 



Tiuj - {E^{N) - En{N)) + iri 

{nN \h{r')\mN) {mN\h{r)\nN) 
nuj + {Em{N)-En{N))+ir] 



(22) 



In the above expression h{r) is the particle density operator defined from the 
wave field operators by 

n{r) = {r)ij{r) (23) 

with ■!/;+(r) = J2k '^k ^k^'"') ^'^d tp{r) = J2k '^fe 0fe('')- By using standard field 
theory techniques it is possible to show that 

X°{r,r';u;;T,)^Yl ^ 'PUr)Mr') G+ir,r' ; Sk + huj;T,) 

k 

+ E /fe MrMir') G;(r, r'; - hu; T,) (24) 

k 

where fk = [1 +exp{(£fc — ^)/fcBTe}]~"^. So far, we have assumed that the 
residual interaction (17) is temperature independent. This assumption is con- 
sistent with the use of a^'j;c('") = c.i;c(") in the calculation of the ground-state 
properties. Therefore, as for = 0, the response function is solution of the 
Dyson equation (16) with given by Eq. (24). 

The above formalism (;an be employed to compute the photo-absorption 
by a metallic nanoparticle of size R. If the wavelength A of the incoming light 
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is such that R the dipolar approximation is valid. From the frequency- 
dependent dipole polarizability 



a (w; Te) = J dn(r; lo; T,) 5Ve^t{r- u) dr (25) 
one obtains the dipolar absorption cross-section [62] 

a(a;;Te) = — Im[a(w;re)]. (26) 

£0C 

As for the zero-temperature case, the dipolar absorption cross-section fulfils 
the well known Thomas-Reiche-Kuhn (TRK) sum rule 



/ 



a (w; Te) du) = . (27) 



Application to femtosecond spectroscopy 

Ultrafast spectroscopy using femtosecond laser pulses is a well suited technique 
to study the electronic energy relaxation mechanisms in metallic nanoparticles 
(see Refs. [6, 12] and references therein). The experiments have been carried 
out with nanoparticles of noble metals containing several thousand atoms and 
embedded in a transparent matrix. By using a time resolved pump-probe con- 
figuration it is possible to have access to the spectral and temporal dependence 
of the differential transmission ^^(r, w), defined as the normalized difference 
between the probe pulse with and without the pump pulse. This quantity 
contains the information on the electron dynamics, which is measured as a 
function of the pump-probe time delay r and of the laser frequency ui. 

For pump-probe delays longer than a few hundred femtoseconds, the ther- 
malization of the electrons is achieved, thus leading to an increase of the 
electron tcmpcratiirc of several hundred degrees. However, the electronic dis- 
tribution is not in thermal equilibrium with the lattice, the thermal relaxation 
to the lattice being achieved in a few picoseconds via electron-phonon scat- 
tering. The energy exchange between the electrons and the lattice can be 
described by the two temperature model leading to a time-dependent electron 
temperature Tf,(t) [63] 

Q^=G{n-Ti), (28) 

where P{t) represents the laser source term, Cj (Ce) is the lattice (electron) 
heat capacity, and G is the electron-lattice coupling factor. In this simplified 
model, the two temperatures arc assumed to be spatially uniform and therefore 
the heat propagation is neglected. 
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Provided that the relative changes of the dielectric function with respect 

to a non-pcrturbed system are weak (linear regime) and that they are only 
due to a modification of the electron temperature, one may identify the spec- 
tral dependence of the differential transmission measured for a given time 
delay as the difference of the linear absorption cross-sections evaluated at dif- 
ferent electron temperatures. More precisely, the differential transmission is 
expressed as 

T ' ' T[Te(0),w] ^ ' ^ ' 

[C7(a;;re(0))-C7(a;;re(r))] (30) 



where / = 2R is the sample thickness (here, the diameter of the nanoparticle) , 
T[Te(r),a;] and T[Te(0),a;] are the probe transmissions in the presence and 
absence of the pump, respectively, and Aa is the pump-induced absorption 
change. Obviously r[T'e(0), w] corresponds to an absorption at room temper- 
ature Te(0) = 300 K for the conditions where the pump-probe experiments 
have been performed. We have computed the optical spectrum of a closed-shell 




2.4 2.8 3.2 3.6 
photon energy (eV) 

Fig. 2. TDLDA photoabsorption cross-section (in atomic units) of Ag2998 encapsu- 
lated in a transparent matrix {em = 1-5) as a function of the photon energy. Solid 
line: Te = K; dashed line: Te = 300 K; dotted line: Te = 1200 K. 



nanoparticle Ag2998 embedded in a transparent matrix (alumina Cm = 1-5) for 
three values of the temperature. The diameter of the nanoparticle is 4.6 nm 
and the photon energy ranges from 2.2 eV to the interband threshold energy 
at 3.8 eV, i.e. in the spectral region associated to the surface plasmon of Ag 

nanoparticlcs. All these values correspond to typical experimental conditions 
performed in our group [6]. The results are presented in Fig. 2. The calculated 
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oscillator strength is 90%. Indeed, due to the presence of the surface plasmon 
resonance, almost all the oscillator strength is concentrated in this energy 
range. A clear red-shift and broadening of the resonance as a function of the 
electron temmperature is observed. 

In the left panel of Fig. 3, the predictions of the normalized differential 
transmission [Eq. (30)] are presented as a function of the photon energy of the 
probe. The comparison is made for two electron temperatures Tg = 600 K and 
Te = 1200 K. The asymmetric shape of AT /T around the resonance energy 
is related to a combination of a red-shift and a broadening of the surface 
plasmon resonance. In the right panel of Fig. 3 the experimental spectrum of 
the normalized AT /T obtained for a pump-probe delay of t = 2ps is depicted. 
The pump pulse is set at 400 nm (second harmonic of a titanium sapphire 
laser amplified at 5 kHz) and the probe comes from a continuum generated in 
a sapphire cristal with the fundamental frequency of the amplified laser [6] . 

The asymmetric spectral shape of the differential transmission spectrum in 
Fig. 3, which is related to the shift and broadening of the plasmon, may have 
several origins. As pointed out in Rcfs. [6, 12, 64], the interband transition 
induces a modification of the real part of the dielectric function in this spectral 
region, the resonance being far enough from the interband threshold to induce 
significant changes of the corresponding imaginary part. As stressed in Refs. 
[12, 64], this is a strong indication that intraband processes also play an 
important role. Indeed, as clearly seen in Fig. 2, the conduction electrons 
contribution leads both to a shift and to a broadening. We can therefore 
conclude that one needs to consider both the interband and intraband part 
on the same footing. Whereas this effect was previously taken into account in 
a phenomenological way via a shifted and broadened Lorentzian shape, here 
we have derived it directly from a quantum many-body approach based on 
the TDLDA at finite temperature. 

2.3 Phase space methods: from Hartree to Wigner and Vlasov 

As we have seen in Sec. 2.1, the most fundamental model for the quantum 
iV-body problem is the Schrodinger equation for the A^-particle wave function 
!Z'(ri, r2, . . . , rjv, i). Unfortunately, the full Schrodinger equation cannot be 
solved exactly except for very small systems. A drastic, but useful and to 
some extent plausible, simplification can be achieved by neglecting two-body 
(and higher order) correlations. This amounts to assume that the A''-body 
wave function can be factored into the product of N one-body functions: 

tZ'(ri,r2, . . . ,riv,i) = il}\{ri,t) ip2{r2,t) . . .'4)N{rN,t). (31) 

For fermions, a weak form of the exclusion principle is satisfied if none of the 
wave functions on the right-hand side of Eq. (31) are identical ^. 

^ A stronger version of the exclusion principle requires that ^(ri , r2, . . . , r]v, t) 
is antisymmetric, i.e. that it changes sign when two of its arguments are inter- 
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photon energy (eV) 

Fig. 3. Left panel; theoretical predictions of the normalized differential transmission 
for Ag2998 embedded in a transparent matrix as a function of the photon energy of 
the probe. Solid line: Te = 600 K; dotted line: Te = 1200 K. Right panel: Normalized 
experimental spectrum of AT/T of silver nanoparticles encapsulated in an alumina 
matrix for a pump- probe delay of 2 ps [6] . 



When the above assumption is made, the A''-body Schrodinger equation 
reduces to a set of one-particle equations, coupled through Poisson's equation 
{time- dependent Hartree model): 

^ ~ ^"^^^ ' = 1 • • • -^orb (32) 

^^=1 (Y.Pc^l^c.f - Mr)] , (33) 



where A'orb > A'' is the number of occupied orbitals, e and m arc the absolute 
electron charge and mass, and e is the dielectric constant; ni{r) is the ion den- 
sity, which is supposed to be fixed and a continuous function of the position 
coordinate. This is known as the 'jcUium' hypothesis and is valid whenever 
the relevant length scales are significantly larger that the ionic lattice spacing 
osiatt ~ 5-^- A-S mentioned in Sec. 2, this is the case for semiconductor nanos- 
tructures, but not so for metals (see Table 1); nevertheless, the jcUium models 
still yields reasonably accurate results for all but the smallest nano-objects. 

The occupation probabilities Pa {J2a=i ~ ■'■) defined to describe a 
Fermi-Dirac distribution at finite electron temperature, Pa = [i- + exp(/3(eQ, — 
IJ,))]~^, where /3 = l/fesTg, /i is the chemical potential, and Eq, is the single- 
particle energy level. In practice, one first needs to obtain the ground-state 
equilibrium solution of Eqs. (32)-(33), which amounts to determining the A^'orb 



changed. This can be achieved by taking, instead of the single product of A'^ wave 
functions as in Eq. (31), a linear combinations of all products obtained by per- 
mutations of the arguments, with weights ±1 (Slater determinant) [65]. This is 
at the basis of Fock's generalization of the Hartree model 
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occupation probabilities and the corresponding energy levels and wave func- 
tions. Subsequently, the equilibrium can be perturbed to study the electron 
dynamics. The numerical methods for the dynamics are quite standard, as the 
Eqs. (32) are basically one-particle Schrodinger equations. We will not enter 
into the details of the numerical methods in this paper: a list of relevant works 
on the Schrodinger equation can be found in Ref. [66]. 

We now show that the Hartree equations can be written in a completely 
equivalent form by making use of the Wigner transformation. The Wigncr rep- 
resentation [29] is a useful tool to express quantum mechanics in a phase space 
formalism (for reviews see [67, 68, 69]). The Wigner function is a function of 
the phase space variables (x, v) and time, which, in terms of the single-particle 
wave functions, reads as 



f{x, V, *) = E ^ P« £^ C + ^, i) - ^, i) e^™"^/'' dX (34) 

(we restrict our discussion to one-dimensional cases, but all results can easily 
be generalized to three dimensions). It must be stressed that the Wigner func- 
tion, although it possesses many useful properties, is not a true probability 
density, as it can take negative values. However, it can be used to compute av- 
erages just like in classical statistical mechanics. For example, the expectation 
value of a generic quantity A{x, v) is defined as: 

/^\ _ J I .f{x,v)A{x,v)dxdv 

~ JJf{x,v)dxdv ' ^"^^^ 

and yields the correct quantum-mechanical value ^. In addition, the Wigner 
function reproduces the correct quantum-mechanical marginal distributions, 
such as the spatial density: 

/-f-00 -^orb 
f{x,V,t) dv =^ Pa \tpa\'^ ■ (36) 
-°° a=l 

We also point out that, of course, not all functions of the phase space 
variables are genuine Wigner functions, as they cannot necessarily be written 
in the form of Eq. (34). In general, although it is trivial to find the Wigner 
function given the wave functions that define the quantum mixture, the in- 
verse operation is not generally feasible. Indeed, there are no simple rules to 
establish whether a given function of x and u is a genuine Wigner function. 
For a more detailed discussion on this issue, and some practical recipes to 
construct genuine Wigner functions, see [70]. 

® For variables whose corresponding quantum operators do not commute (such as 
xv), Eq. (35) must be supplemented by an ordering rule, known as Weyl's rule 
[69]. 
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The Wigner function obeys the following evolution equation: 

dt dx 



A\ . / A 



-n""2 



(37) 

where (^(x, t) is the self-consistent electrostatic potential obtained self-consistently 

from Poisson's equation (33). 

Developing the integral term in Eq. (37) up to order Oih?) we obtain 

^+^,^ + ^^^^^^^ + 0(^4). (38) 

dt dx m dx dv 24,m^ dx^ dv^ 

In the limit ?l ^ one recovers the classical Vlasov equation, well-known from 
plasma physics (see Fig. 1). The Vlasov-Poisson system has been used to study 
the dynamics of electrons in metal clusters and thin metal films [45, 18]. It is 
appropriate for large excitation energies, for which the electrons' de Broglie 
wavelength is relatively small, thus reducing the importance of quantum ef- 
fects in the electron dynamics. Nevertheless, for metallic nanostructures at 
room temperature, the equilibrium must be given by a Fermi-Dirac distri- 
bution, because the Fermi temperature is very high (see Table 1). For semi- 
conductor nanostructures, Tp ^ 10 — 50K, so that a Maxwell-Boltzmann 
equilibrium is sometimes appropriate for moderate temperatures. 

The Wigner equation must be coupled to the Poisson's equation for the 
electric potential: 

-^ = -^[ni{x)-n{x,t)], (39) 

The resulting Wigner-Poisson (WP) system has been extensively used in the 
study of quantum transport [71, 72, 73]. Exact analytical results can be ob- 
tained by linearizing Eqs. (37) and (39) around a spatially homogeneous equi- 
librium given by no/o(f) (Maxwell-Boltzmann or Fermi-Dirac distribution), 
where no = rij = const, is the uniform equilibrium density. By expressing the 
fluctuating quantities as a sum of plane waves exp(iA;a; — icot) with frequency 
u) and wave number k, the dispersion relation can be written in the form 
e{k, (j) = 0, where the 'dielectric constant' e reads, for the WP system, 

.wp(-, k) = i + —^j nn^j^) 



or equivalently 



fo{v) 



ewp(u;,/c) = 1 - / , ' \2.4M 2 . (41) 

no J [u — kvy — n fc*/4rn^ 
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This is just the Lindhard [74] dispersion relation, well known from solid-state 
physics. From Eq. (40), one can recover the Vlasov-Poisson dispersion relation 
by taking the classical limit 

eM-,k) = l + ^J^dv. (42) 

nok J uj — kv 

The equivalence of the Hartree and Wigner-Poisson methods can be easily 
proven by comparing the linear results. For the Hartree equations (32), we 
linearize around a homogeneous equilibrium given by plane waves: 




each with occupation number p„ and energy = muQ^/2. The Hartree di- 
electric constant is found to be 

en{u;,k) = 1 - J^p„ , ^f" ,,,, (44) 

which is a discrete form of the Wigner-Poisson dispersion relation (41). 

Example — Ultrafast electron dynamics in thin metal films 

Several experiments have shown [2, 3] that electron transport in thin metal 
films occurs on a femtosecond time scale and involves ballistic electrons trav- 
eling at the Fermi velocity of the metal vp- More recently, a regime of low- 
frequency nonlinear oscillations (corresponding to ballistic electrons bouncing 
back and forth on the film surfaces) was measured in transient reflection ex- 
periments on thin gold films [75]. 

These findings were corroborated by accurate numerical simulations based 
on the one-dimensional Vlasov-Poisson equations [18]. The electrons are ini- 
tially prepared in a Fermi-Dirac equilibrium at finite (but small) tempera- 
ture. They are subsequently excited by imposing a constant velocity shift 
Av = O.OSvp to the initial distribution, which is a rather strong excitation. 
This scenario is appropriate when no linear momentum is transferred parallel 
to the plane of the surface (i.e., q\\ = 0) and is relevant to the excitation of 
the film with optical pulses [76]. For = 0, only longitudinal modes (volume 
plasmon with u> = Up) can be excited. 

As a reference case, we studied a sodium film with initial temperature Tg = 
O.OOSTp ~ 300 K and thickness L ~ 120 A. The time evolution of the thermal 
i?th and center-of-mass i^cm energies was analyzed (Fig. 4). During an initial 
rapidly-oscillating phase, £^cm is almost entirely converted into thermal energy 
(Landau damping). After saturation, a slowly oscillating regime appears, with 
period equal to 50uj~^ ~ 5.3fs, where LOp = {e^n/meoy^^ is the plasmon 
frequency. This period is close to the time of flight of electrons traveling at 
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Fig. 4. Time evolution of the thermal, potential, and center-of-mass energies of the 
electron population in a thin sodium film. 



the Fermi velocity and bouncing back and forth on the film surfaces (further 
details are provided in our previous work [18]). 

The phase space portrait of the electron distribution, shown in Fig. 5 
clearly reveals that the perturbation starts at the film surfaces, and then 
proceeds inward at the Fermi velocity of the metal. The structure formation 
at the Fermi surface, which has spread over the entire film for LOpt > 150, 
is responsible for the increase of the thermal energy (and thus the electron 
temperature) observed in Fig. 4. As no coupling to an external environment 
(e.g., phonons) is present, this excess temperature cannot be dissipated. 
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Fig. 5. Phase space portrait of the electron distribution. Velocity is normalized to 
the Fermi velocity, and space to the Thomas- Fermi screening length Lf = vp/ujp. 
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Quantum simulations of the electron dynamics using the Wigner-Poisson 
system were performed more recently: as expected, the Vlasov results were 
recovered in the large excitation regime 6v > Q.OSvp. For smaller excitations, 
a different regime appears, in which the ballistic oscillations described above 
are no longer observed. Further work is in progress on this issue [77] . 

Beyond the mean field 

The mean-field approach described above is appropriate to describe the elec- 
tron dynamics on very short time scales (< f OOfs). On a longer time scale (O.f- 
Ips), the injected energy is redistributed among the electrons via electron- 
electron (c-e) collisions. Elcctron-phonon (e-ph) thermalization (i.e., coupling 
to the ionic lattice) is generally supposed to occur on even longer time scales. 
However, the results of Refs. [5, 78] on thin gold films have shown that 
nonequilibrium electrons start interacting with the lattice earlier than ex- 
pected, so that a clear-cut separation between e-c and c-ph relaxation is not 
entirely pertinent. 

The phase-space approach is particularly well-suited to include corrections 
that go beyond the mean-field picture. This can be done with relative ease for 
semiclassical models (Vlasov) , by using a Boltzmann-like e-e collision integral 
that respects Pauli's exclusion principle (Uhling-Uhlenbeck model) [79] : 

(I),, = / W -(^)l^-l(^^/^7374 - hfJJ.) , (45) 

where Vi2 is the relative velocity of the colliding particles 1 and 2, a{f2) is the 
differential cross section depending on the scattering angle f2, and indices 3 
and 4 label the outgoing momenta, fi = f{r,pi,t) and /j = 1 — /i/2. This 
collision term is similar to the well known classical Boltzmann collision term 
but for Pauli blocking factors fifj- As known from solid-state physics, this 
blocking factor plays a dramatic role for electronic systems [65]. At Te = 
OK, all collisions are Pauli blocked and the coUisional mean-free path of the 
electrons becomes infinite. But if the system becomes excited, phase space 
opens up and activates the collision term. The effect of the above e-e collision 
term on the semiclassical Vlasov dynamics in metal clusters was investigated 
numerically in [80]. 

It is conceptually harder to include collisions in fully quantum models. A 
significant constraint is that nonunitary corrections to the Wigncr equation 
should be written in 'Lindblad form' [81], which guarantees that the evolved 
Wigner function corresponds to a positive-definite density matrix. 

The Uhling-Uhlenbeck collision term (45) is a complicated nonlinear inte- 
gral, which is difficult to implement in a numerical code. It is therefore useful 
to construct some simplified collision terms that are more easily amenable to 
numerical treatment. In the following, we briefly illustrate two simple models 
of e-e and e-ph collisions that we have employed in our previous works. 
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Electron- electron collisions. 

To model e-e collisions, a relaxation term is added to the right-hand side of 
the Vlasov or Wigner equation: 

= -l/ee(Te)(/-/oo), (46) 

e— e 

where i/ge is the average e-e collision rate and f^ (x,v) is a Fermi-Dirac dis- 
tribution. The idea behind this model is that the electron distribution will 
eventually relax, on a time scale of the order towards a Fermi-Dirac 
equilibrium with total energy equal to that of the initial electron distri- 
bution f{x,v,t = 0+), including of course the initial excitation energy. For 
electrons near the Fermi surface, the e-e collision rate can be written as [82] 

^/ee(7;) = a(fcBTe)^ (47) 

where a is a (dimensional) proportionality constant. The latter has been esti- 
mated from numerical simulations of the electron dynamics in sodium clusters 
[80], yielding a ~ 0.4 fs~^eV~^, which is also compatible with the analytical 
prediction given by the random phase approximation [82] . The electron tem- 
perature is computed instantaneously during the simulation, and plugged into 
the expression for the collision rate (47). It is important to underline that the 
above model for e-e collisions, though simple, is completely self-contained and 
requires no additional ad-hoc parameters. The model has been applied to the 
electron dynamics in thin metal films. The slow ballistic oscillations of Fig. 4 
are still observed, although they are damped on a time scale of the order of 
SOOw-gi ~ 50fs (see Fig. 6). 
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Fig. 6. Evolution of the thermal energy for a case with e-e collisions and L = 
IOOLf ^ 120A. 



Electron-phonon collisions. 

By coupling to the ionic lattice, the electrons progressively relax to a thermal 
distribution with a temperature equal to that of the lattice Tj. This relaxation 
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time is generally termed ti in the semiconductor literature. In addition, the 

lattice acts as an external environment for the electrons, leading to a loss of 
quantum coherence over a time scale (decoherence time). The relaxation 
and decoherence times correspond, respectively, to the decay of diagonal and 
nondiagonal terms in the density matrix describing the electron population. 

Such environment-induced decoherence can be modeled, in the Wigner 
representation, by an appropriate friction-diffusion term [83] : 



where 7 is the relaxation rate (inverse of the relaxation time n), and Dy, 
Dx are diffusion coefficients in velocity and real space respectively, which are 
related to the decoherence time T2 and depend on the lattice temperature Tj. 
The effect of the diffusive terms is to smooth out the fine structure of the 
Wigner function, thus suppressing interference phenomena, which are a typi- 
cally quantum effect. Finally, we recall that, in order to preserve the positivity 
of the density matrix associated to the Wigner distribution function, the e- 
ph collision term (48) must be in Lindblad form [81]. This is automatically 
achieved [84] if the coefficients respect the inequality D^Dx > ^"^Ti^ j^rr?. 

2.4 Hydro dynamical models: from micro to macro 

Despite its considerable interest, the Wigner-Poisson (WP) formulation presents 
some intrinsic drawbacks : (i) it is a nonlocal, integro-differential system; and 
(ii) its numerical treatment requires the meshing of the whole phase space. 
Moreover, as is often the case with kinetic models, the Wigner-Poisson sys- 
tem gives more information than one is really interested in. For these reasons, 
it would be useful to obtain an accurate reduced model which, though not 
providing the same detailed information, is still able to reproduce the main 
features of the physical system under consideration. 

In this section, we will derive an effective Schrodinger-Poisson (SP) sys- 
tem, which, in an appropriate limit, reproduces the results of the kinetic WP 
formulation [85] . In order to obtain the effective SP system, we will first derive 
a system of reduced hydrodynamic (or fiuid) equations by taking moments of 
the WP system. It will be shown that the pressure term appearing in the fluid 
equations can be decomposed into a classical and a quantum part. With some 
reasonable hypotheses on the pressure term, the fluid system can be closed. 
For simplicity of notation, only one-dimensional problems will be considered, 
but the results can be easily extended to higher dimensions. 

In order to derive a fiuid model, we take moments of Eq. (37) by integrat- 
ing over velocity space. Introducing the standard definitions of density, mean 
velocity, and pressure 




(48) 



n= f dv , u = — fv dv , P = m 





(49) 
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it is obtained 

du du e dS 1 dP 

^+^^ = -^ ^1 

at ox max mn a x 

We immediately notice that, surprisingly, Eqs. (50)-(51) do not differ from the 

ordinary evolution equations for a classical fluid. It can be shown, however, 
that quantum effects are actually hidden in the pressure term, which may be 
decomposed into a classical and a quantum part. 

By using the definition of the Wigner function (34) and representing each 
state in terms of its amplitude ^Jn^ and phase Sa 



ipa{x,t) = y/na{x,t) exp {iSa{x,t)/h), (52) 
we obtain that P = P'-^+P'^. The classical part of the pressure can be written 



as 



E lla 2 I \ " 
Pa - l^Pa U. 



= mn{{ul) - {u^f), (53) 




where mua = dSa/dx [the u^s shoiild not be mistaken with the global mean 
velocity u defined in Eq. (49)] . This is the standard expression for the pressure 
as velocity dispersion, thus justifying the term 'classical' pressure. 
The quantum part of the pressure is written as 



(54) 

It can be shown that, for distances larger that the Thomas- Fermi screening 

length Lp, one can replace with n, the total density as defined in Eq. (49). 
In order to close the fluid system (50)- (51) one still has to express the classical 
pressure in terms of the density n. This is the standard procedure adopted 
in classical hydrodynamics: the relation P'^{n) is the equation of state, and 
depends on the particular conditions of the system, notably its temperature. 
With these hypotheses, the Eq. (51) reduces to 

du ^ ^du _ e 8(1) 1 dW ^ d / d^{y^)/dx^ \ 
dt dx mdx mdx 2m^ dx \ \fn J ' 

where we have defined the effective potential 

Equations (50) and (55) constitute the quantum hydrodynamical approxima- 
tion to the full Wigner (or Hartree) equation. 
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It is now possible to combine Eqs. (50) and (55) into an effective nonlinear 
Schrodinger equation. To this purpose, let us define the effective wavefunction 



= ^n{x, t) exp {iS{x, 1)/%) , (57) 

with S{x,t) defined according to mu{x,t) = dS{x,t)/dx. We obtain that 
^{x, t) satisfies the equation 

By linearizing Eqs. (50) and (55) around a homogeneous equilibrium, we 
obtain the following dispersion relation 

a;^=a;^+Vfc^ + |^, (59) 

where tuVq = {dP'~^ /dn)n=no- It can be proven that, by an appropriate choice 
of the equation of state P^{n), Eq. (59) reproduces correctly the leading terms 
of the Hartree or Wigner dispersion relation. 

To summarize, we have shown that, under appropriate conditions, the 
Hartree or Wigner models can be reduced to a set of two hydrodynamical 
equations (50) and (55), or, equivalcntly, to a single nonlinear Schrodinger 
equation (58). The two hypotheses used for this reduction were that (i) all 
quantities vary on a length scale larger than Lp; and (ii) the equation of state 
for the classical pressure is P*^ = P^{n) (standard fluid closure). 

Example — Thin metal films 

We have studied the electron dynamics in a thin metal film using the above 
quantum hydrodynamical model [86]. A preliminary result is shown in Fig. 
7, where we plot the evolution of the thermal and potential energies against 
time. In order to compare to the Vlasov simulations described in Sec. 2.3, 
the hydrodynamic equations arc solved in the scmiclassical limit, i.e. using 
a small value of the Planck constant normalized to Ep/ujp (note however 
that here the initial excitation 5v = 0.22vf is larger compared to the case of 
Fig. 4, where Sv = 0.08?;^?). The hydrodynamic results display some coherent 
oscillations at high frequency, which are a typical signature of quantum effects. 
Nevertheless, the initial increase of the thermal energy is clearly captured 
and the subsequent ballistic oscillations are still visible, particularly on the 
potential energy. 



3 Spin dynamics 



The dynamics of magneto-optical processes in metallic nanostructurcs de- 
pends on the temporal and spatial characteristics that are being investigated. 
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Fig. 7. Time evolution of the thermal and potential energies (normalized to Ef) of 
the electron population, obtained using a quantum hydrodynamics model. 



Short time scale (t < 10^^^ s) have only been explored recently. In 1996, 
the group of Jean- Yves Bigot in Strasbourg highlighted the existence of ul- 
trafast demagnetization processes (within less than a hundred femtoseconds) 
induced by femtosecond laser pulses in ferromagnetic thin films [24]. These 
demagnetization processes are not yet fully understood. 

Prom a theoretical point of view, very little is known on the time-dependent 
magneto-optical response of metallic nanostructures to an ultrafast optical 
pulse. The main difficulty is to provide an adequate description of the interplay 
between electronic and spin degrees of freedom in the metal. So far, only two 
theoretical models have been proposed to explain this effect [25, 26]. These 
works are based on two different mechanisms: in [25], the spin-orbit coupling 
is invoked, whereas in [26] phonon or impurity mediated spin-flip scattering 
is privileged. Unfortunately, the parameters employed in [25] are not realistic 
and the model developed in [26] is a phenomenological approach that does 
not allow quantitative predictions. Prom the above considerations it follows 
that there is a need for the development of efficient theoretical models able to 
explain in a quantitative manner the experimental findings. 

A proper treatment of spin dynamics requires an extension of our model 
(TDLDA) to include; spin degrees of freedom. In the following, the formalism of 
the time-dependent local-spin-density approximation (TDLSDA) in the linear 
regime (including also its extension to finite temperature) is presented. A 
second part will be devoted to the nonlinear dynamics. 

3.1 Linear response: local spin density approximation 

The generalization of the linear TDLDA to spin polarized electron systems 
has been performed by Rajagopal [87]. In the following we provide the basic 
equations of this approach including its extension to finite temperature. 

Within the framework of DPT one can calculate the spin density matrix 
fiacr' (f) defined as 
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n,,,{r) = {0\i;+ir)i;,,ir)\0) (60) 

where tp'^{r) and V'ct(^) are the wave field operators corresponding to the 
creation and annihilation of an electron with spin a at position r, and |0) is the 
ground state of the system. When the system is subjected to a small local spin- 
dependent external potential 6V^Jl {r;uj) (this quantity describes the coupling 
of the charge and spin of the electrons to external electric and magnetic fields) 
the spin-density response function is defined through the equation 

Sn,,, {r;u;)=J2 / Xaa',.,a, {r, r'; oj) J^f^ (r'; u;) dr' . (61) 

(Tia'2 

For the sake of simplicity, we restrict ourself to the case of coUinear mag- 
netism, i.e. to the case of a uniform direction of magnetization. This restric- 
tion leads to a diagonal spin-density matrix (n^a' = n^S^cr') and simplified 
expressions. The spin-density response function defined in Eq. (61) reduces to 

6n,{r;u;) = ^ / Xa.'{r,r';oj) SV^^^ir'-uj) dr' (62) 

which can be rewritten as 

5n,(r; co) = J2 [ xL' (r, r'; to) to) dr' (63) 

with 



e^/47reo 



--^+f^f{r,r';u;)\6n^.{r';u;)dr'. (64) 



In the above expression the function f^f {r,r';oj) is the Fourier transform of 
the time-dependent kernel defined by /^^ {r,t;r',t') = 5V^^{r,t)/Sna-i{r',t') 
and x°^/(r,r';a;) is the non- interacting retarded spin-density correlation 
function. For spin polarized electron systems the exchange-correlation po- 
tential is defined as 



d 

dUa 



{nwa:c(n+,n_)} 



(65) 



n_|_=n_|_(r);n_=n_ (r) 



where f2xc[n+,ri^] = J n(r)ujxcin+{r),n^{r)) dr is the exchange-correlation 
thermodynamic potential and u>xc the exchange-correlation thermodynamic 
potential per particle of the homogeneous electron gas calculated at the local 
density n and magnetization m = n+ — n_ . By noting that 



d 



{nwa:c{n+,n-)} 



d_ 
dn 



{nujxc{n,m)} 



{nujxc{n,m)} , 
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the expression (65) can be rewritten as [88] 



+ anBBa:c{r) , (66) 

n=n{r);m=m{r) 

where B^c{r) = Ms^ {ncJa^dn, m)}]^^^^^y^^^_^^^^ is the exchange-correlation 
magnetic field acting on spin, and hb = eh/ (2m) is the Bohr magneton. This 
is an internal magnetic field. The response functions x° and x are related by 
an integral equation (to be more precise, due to the spin degree of freedom, 
it is a matrix integral equation) 

Xaa' {r, r'; w) = xL' (r, r'; w) + ^ J J xla, (r, r"; u) 

(71(72 

X K-'--{r",r"';u) Xa,a'ir"',r';u;) dr"dr"', (67) 
with the residual interaction defined by 

if--(r,r';a;) = _-£!_^<5,^,, + /--(r.,r>). (68) 

As for TDLDA, in the adiabatic local-density approximation (ALDA) the 
exchange-correlation kernel is frequency-independent and local and reduces 
to 

6{r-r') . (69) 

n=n{r);m=Tn{r) 

It should be mentioned that the functional w^c in the above expression should 

be the same as the one used in the calculation of the groimd state (sec 
Eq. (65)). By using the same field-theory techniques employed previously for 
TDLDA (see Sec. 2.2), one can show that the free response function reads 

x'i,,{r,r';Lo;T,)=S,,,Yl 4>l*{r)ct>l{r') Gl{r,r';el + hLo;Te) 

k 

+ J2fk rk{rWk*{r') Gl*{ry;e% - nu:;^) , (70) 
fe 

where 4)%{r) and ef, are the one-electron Kohn-Sham wave functions and en- 
ergies, respectively. GZ^ is the one-particle retarded Green's function for the 
spins (J and = [l + exp {{el - iJi)/kBTe)]~ . Similarly to TDLDA, we have 
assumed that the residual interaction (68) is temperature independent. Thus, 
it is consistent with the use of u:xc{n, m) = e^dn, m) in the calculation of the 
ground-state properties. 

From the above formalism one can compute the dipolar absorption cross- 
section 

a(a;;Te) = — Im[a(a;;Te)], (71) 



— {nuja:c{n,m)} 



J xc 



\r,r') 



c(n,TO)] 



Collective Electron Dynamics 31 



where a is the frequency-dependent dipole electric polarizability defined as 

a{uj;T^) = j [5n+{r;Lo;T^)+5n_{r;uj;T^)] SV^^tir;'^) dr . (72) 

By analogy, one defines a quantity which is constructed from the local mag- 
netization (instead of the local density) 

am (w; Te) = — Im [am {to; Te)] , (73) 

SqC 

where am is the frequency-dependent dipole magnetic polarizability defined 



as 



am {to;Te) = j [dn+{r;u;Te) - 5n-{r;Lo;Te)] SVe^tir;u) dr . (74) 
On can show that am fulfils the following sum rule 

am (u;;r,)d^= ^^'^^^-^ (75) 
c 

where M = — N~ is the total magnetization of the system (A'^^ being the 
number of spins up and N~ the number of spins down). It is worth mentioning 
that M is generally temperature dependent [89]. 

3.2 Nonlinear response: Phase-space methods 

In order to investigate the nonlinear regime of the charge and spin dynamics, 
a phase-space approach is particularly interesting. In this paragraph, wo will 
construct a Wigncr equation that includes spin effects in the local density 
approximation, and show that its classical limit takes the form of a Vlasov 
equation. 

The starting point for the derivation are the time-dependent Kohn-Sham 
(KS) equations described in Sec. 3.1. In terms of the Pauli 2-spinors 



!^,'(r,t) 



'^^/(r,t) 



the KS equations can be written as: 

W,{r,t) (76) 



-^V^ + V{r,t) ]I + i,B<y ■ B{r,t) 



where V{r,t) = Vcxi{r,t) + VH{r,t) + V^" (r,t), is Bohr's magneton, cr = 
(cx) fj/) fz) are the 2x2 Pauli matrices, and I is the identity matrix. Here, 
Kixt is an external potential (e.g. ionic jellium, external electric field, ...), Vh 
is the Hartreo potential that obeys Poisson's equation, and Kf^, is the scalar 
part of the exchange-correlation potential. The magnetic field B = Bext + -Bxc 
is composed of an external part and an 'internal' part that stems from the 
exchange and correlation energy [see Eq. (66)]. In the so-called 'coUinear' 
approximation, the latter reduces to Bxc = -Bxc-^- 
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Equation of motion for the density matrix 

By defining the density matrix 

p^^'{v,v')=Y,^^{r)^f*ir') (77) 

i 

where rj =t, i, the KS equations (76) can be written in the fohowing compact 
form (Von Neumann equation): 

ih^ = [H,p], (78) 

where 

The only nondiagonal terms in the Hamiltonian come from the external or 

internal magnetic field B. 

We now introduce the following basis transformation for the Hamiltonian: 

H = hoI + h-a (80) 

where h = {h^, hy,hz), and 

ho = ^ , = ^ (81) 

hz = ^ , hy = (82) 

For the Hamiltonian of Eq. (76) , we have 

ho{r) = ^^W^+V{r,t) (83) 
ha{r) = I^BBa{r,t), a = x,y,z (84) 

The same transformation (with identical notation) is also applied to the den- 
sity matrix. With these definitions, the equations of motion for po and pa read 
as 

ihdtpo = [ho,po]+ ^ [ha, pa] (85) 

a—x,y,z 

ihdtPa = [ho, Pa] + [ha, Po]- (86) 



'Spin' Wigner and Vlasov equations 

By making use of the Wigner transformation 
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A 



Tfl f 

one can easily obtain the equations of motion for the Wigner functions: 



(87) 
(88) 



^/o+v-/o- 



m 



4^ 2iTrn^ J J 



im(v—v'^X/h 



/o(r,v',t)- 

/„(r,v',i)=0 



[dX [ dv'e^'"(^--')^/^ 
2iTrn^ J J 

^ [dX [ dv'e'-(---')V^ 
2inh^ J J 



Br, 



X 

^+2 



/a(r,v',t)- 

/o(r,v',t) = 



The corresponding Vlasov equations are obtained in the classical limit ?l — > 0: 

d_ 



-Jo + V— /o 7^ — 7^ > — — = 



9r " 



9r d\ 



dr dv 



If +v— / -l^^-^^^=0 
9t " 9r " m dr dv m dr dv 



(90) 



with a = x,y, z. 

Within the collinear approximation, the equations for a = x,y vanish. In 
this case, it is more convenient revert to the original representation and use 

/t = /o + 

/i = /O - fz- 

The corresponding Vlasov equations then become 



ot or m \ or or J ov 

If +v— / 1 (dV i^j^^^A -0 
dt dr ml dr dr J dv 



(91) 
(92) 



The above Wigner and Vlasov equations can be used to study the nonlin- 
ear spin dynamics in a ferromagnetic nanoparticlc or thin film, using numer- 
ical techniques similar to those employed for the electron dynamics. In their 
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present form, these equations preserve the total spin, and thus cannot be used 
to describe the loss of magnetization observed in experiments [24]. A proper 
generalization, along the lines of the e-e and e-ph collision operators detailed 
in Sec. 2.3, would be necessary to account for these effects. 

4 Numerical example: the nonlinear many-electron 
dynamics in an anhcirmonic quantum well 

In order to illustrate qualitatively the practical implementation of the models 
described in the previous sections, we concentrate on a specific - and rela- 
tively simple - example. We consider an electron population confined in a 
one-dimensional anharmonic well defined by the potential 



where to,k is the effective electron mass. The frequency wq can be related 
to a fictitious homogeneous positive charge of density no via the relation 
u)q = e^no/m^e. The total potential seen by the electrons is the sum of the 
confining potential Konf and the Hartree potential, which obeys Poisson's 
equation 



where e is the absolute electron charge and e is the effective dielectric constant. 
As initial condition, we take a Maxwell-Boltzmann distribution with Gaussian 
density profile 



with temperature and peak density rie- 

The electron dynamics is mainly determined by two dimensionless param- 
eters: (i) the 'filling fraction' r] = rie/no = uip/uiQ, which is a measure of 
self-consistent effects (in the limit case rj = 0, corresponding to very dilute 
electron densities, the Hartree potential is negligible); and (ii) the normalized 
Planck constant H = Jiwo/kBTe, which determines the importance of quan- 
tum effects. Notice that a small value of H corresponds to a large electron 
temperature. 

We use typical parameters for semiconductor quantum wells [90, 91]: ef- 
fective electron mass and dielectric constant m* = 0.067rne and e = 13eo; 
volume density uq = 10"'^^cm~^, oscillator energy hujQ = 3.98mcV, and oscil- 
lator length Lho = \/h/m^LOo ~ 17nm. For ry = 1, this yields a maximum 
surface density for the electrons Ug = 4.64 x 10^°cm~^ and a maximum Fermi 
temperature Tp = 29. 3K. A low electron temperature Tg ~ 46K then yields 
~ 1, whereas at room temperature Tg ~ 300K one has ~ 0.15. 



Konf(a;) = -Lolm^x'^ + -Kx'^ 



(93) 




(94) 




(95) 



Collective Electron Dynamics 35 



The electron dynamics is excited by shifting the electron density of a fi- 
nite distance Sx ~ L\^o- We will primarily be interested in the relaxation 
of the electric dipole, defined as the center of mass of the electron pop- 
ulation: d{t) = J J fxdxdv/ J J fdxdv, and of the average kinetic energy 
■E'kin ^ J J fiTi^v^dxdv/ J J fdxdv. 

First, we present results obtained from the numerical resolution of the 
Wigner equation (37), coupled to Poisson's equation (94). The results were 
obtained with a numerical code that combines the split-operator method with 
fast Fourier transforms in the velocity coordinate [92]. We explore the electron 
dynamics for different values of the two relevant dimensionless parameters, H 
and T]. The anharmonicity parameter appearing in the confining potential (93) 
is fixed to K = 0.1 (in units where U = m* = lvq = 1). If the confinement 
were purely harmonic (i.e., K = 0), the dipole would simply oscillate at the 
frequency luq irrespective of the value of the filling fraction. This result goes 
under the name of Kohn's theorem [31], and we have checked that it holds for 
our numerical simulations. When the confinement is not harmonic, the dipole 
should decay because of phase mixing effects. 

The numerical results are shown in Fig. 8 (dipole) and Fig. 9 (kinetic 
energy). The fast oscillations correspond to the center of mass of the electron 
gas oscillating in the anharmonic well. For low electron densities and large 
temperatures [r] = 0.1, Tg = 300K), the dipole relaxes to the bottom of the 
well, d ~ 0, whereas the kinetic energy relaxes to a constant asymptotic value. 
This is a semiclassical regime where the energy spectrum is almost continuous: 
the observed relaxation is due to phase mixing effects. 

Decreasing the temperature {T^ = 46K) while keeping the density low 
(ji = 0.1) produces a revival that occurs after the kinetic energy has initially 
relaxed. This is a typically quantum effect resulting from the discrete nature 
of the energy spectrum. The revival is clearly visible on the kinetic energy, but 
not so much on the dipole. When the electron density is large (77 = 1), self- 
consistent electron-electron interactions (Hartree potential) prevent the dipole 
and the kinetic energy from relaxing completely, even at large temperatures. 

Next, we have added a dissipative term to the Wigner equation, in order 
to model electron-phonon (e-ph) collisions. This model has been discussed in 
Sec. 2.3. The relaxation rate is chosen to be 7 = O.OOlwo, yielding a realistic 
relaxation time ri = ~ 165ps. The velocity-space diffusion coefficient is 
Dy = jvth, where the thermal velocity is Vth = \fk^T^Jr(i^. The relaxation 
time T2 depends on the velocity scale: for instance, a velocity scale Av is 
damped on a time scale T2 = t\Av jvih. Therefore, for velocity scales smaller 
than the thermal velocity, the decoherence time is always smaller than the 
relaxation time, in accordance with experimental findings. 

We simulated the low temperature scenario (Tg = 46K) in the presence of 
e-ph collisions, and observed that the revival occurring in the kinetic energy 
for 7/ = 0.1 is now suppressed (see Fig. 10). For large densities, however, the 
coherence of the electron motion is not lost, and the relaxation of the dipole 
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Fig. 8. Evolution of the electric dipole (in units of Lho = 17nm) obtained from the 
Wigner-Poisson model, for several values of rj and the electron temperature. Time 
is normalized to the oscillator frequency. 



and the kinetic energy is only marginally faster compared to the collisionless 
regime. 

Finally, we want to consider the zero-temperature case. For doing this, 
we resort to the hydro dynamical model described in Sec. 2.4. The relevant 
dimcnsionlcss parameters now arc rj and rgo, the normalized Wigner-Seitz 
radius computed with the background density uq. For no = lO^^cm"^, one 
has = 2.8. In Fig. 11 we plot the evolution of the electric dipole for 
different values of the filling fraction. Now, even for low electron densities, 
the dipole oscillates indefinitely without any appreciable decay. For larger 
electron densities, the motion is even more regular. It appears, therefore, that 
the dynamics becomes more and more regular as the electron temperature 
decreases, i.e. when quantum effect become more important. As mentioned 
above, this is essentially due to phase mixing effect, which become increasingly 
important in the semiclassical regime, where the energy levels are almost 
continuous. 



5 Conclusions and Perspectives 



In this review paper, we have presented some of the most common theoret- 
ical models used to describe the charge and spin dynamics in metallic and 
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ri=l ; Te = 300K 7? = 0.1 ; Te = 300K 




Fig. 9. Evolution of the kinetic energy (normalized to fiujo = 3.98meV) obtained 
from the Wigncr-Poisson model, for several values of and the electron temperature. 
Time is normalized to the oscillator frequency. 



semiconductor nanostructures. Three levels of description have been identi- 
fied (sec Fig. 1): (i) the full quantum TV-body problem, which can only be 
addressed for small systems by using, for instance, the Configuration Inter- 
action (CI) method; (ii) mean field models (Hartree and Wigner) and their 
generalizations to include exchange and correlations (Hartree-Fock, density 
functional theory); and (iii) quantum hydrodynamical models, which describe 
the electron dynamics via a small number of macroscopic variables, such as 
the density and the average velocity. 

Each of these quantum-mechanical approaches has its classical counter- 
part: classical A/'-body models have been developed for molecular dynamics 
computations, as well as for gravitational iV-body problems; classical mean 
field models are ubiquitous in plasma physics (Vlasov-Maxwell equations) and 
in the study of self-gravitating objects such as star clusters, galaxies, or even 
the entire universe; classical hydrodynamics hardly needs mentioning, as it is 
in itself an extremely wide field of research. 

For each approach, we have stressed the difference between the linear and 
the nonlinear response. The former is valid for weak excitations and pre- 
supposes that the response is directly proportional to the excitation. Linear 
response theory is generally represented in the frequency domain. In contrast, 
nonlinear effects kick in for large excitations, and are best described in the 
time domain (this is because the time-frequency Fourier transform is a linear 
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Fig. 10. Evolution of the kinetic energy (top panels) and electric dipole (bottom 
panels), from the Wigner-Poisson model including e-ph collisions. Same normaliza- 
tions as in Figs. 8 and 9. 
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Fig. 11. Evolution of the electric dipole for r; = 1 (left frame) and q = 0.1 (right 
frame), obtained from the quantum hydrodynamic model at Tg = 0. 
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operation, thus not adapted to describe nonlinear relations). Although a vast 
literature on the linear electronic response is available and dates back from 
the works of Drude in the early twentieth century, nonlinear effects have only 
been investigated in the last two decades, mainly with computer simulations. 

The mean field level of description is perhaps the most widely used, as it 
incorporates, at least to lowest order, some of the features of the A^-body dy- 
namics, but still avoids the formidable complexity of the full problem. A par- 
ticularly challenging open problem is the inclusion of dynamical correlations 
within mean-field models. Dynamical correlations differ from the correlations 
that are included in time-dependent density functional theory (TDDFT), inas- 
much as they cannot be described by a slowly-varying density functional, as 
is done in ALDA (adiabatic local-density approximation). Whereas adiabatic 
correlations are described within an essentially Hamiltonian formulation and 
thus cannot model irreversible effects, dynamical correlations are responsible 
for the relaxation of the electron gas towards thermodynamical equilibrium. 
Some recent results have been obtained using a generalization of TDDFT that 
relies on the electron current as well as the electron density [42]. The phase- 
space approach, via the Wigner formulation, also appears promising to model 
effects beyond the mean field, as we have illustrated in Sec. 2.3. 

Another important issue, which was not mentioned earlier in this review, 
is the inclusion of relativistic corrections in the above models for the electron 
dynamics. Spin-orbit coupling (which is an effect appearing at second order in 
v/c) is sometimes taken into account in a semi-phcnomcnological way within 
the Pauli equation. However, other terms occurring at the same order are often 
neglected without further justification. A consistent derivation of relativistic 
effects to a certain order in v/c can of course be carried out, starting from the 
Dirac equation, for the case of a single particle in an external electromagnetic 
field [93]. For a many-body system, this issue is much trickier and is the object 
of current investigations. 

Nanostructures are by definition finite-size objects. Due to the presence of 
boundaries and interfaces, the electron dynamics can thus display novel and 
unexpected features compared to bulk matter. For example, as the elastic 
and inelastic scattering length (~ 10 — 50nm for bulk metals) are much longer 
than the size of the system, an electron - or a group of electrons - can travel 
coherently through the length of the system, thus leading to ballistic transport 
between the surfaces. The theoretical tools to study finite-size nano-objects 
are also relatively recent, and have been developed alongside the experimental 
breakthroughs that made these objects widely available. 

If the electron dynamics in nanosized objects has received considerable at- 
tention for the last thirty years, the spin dynamics is a much younger field of 
research, both experimentally and theoretically. Nevertheless, the already ex- 
isting applications to memory storage and processing, and the still speculative, 
but highly enthralling, developments in quantum computing, have stimulated 
a large number of works in this direction. In Sec. 3 we have illustrated how the 
models for the electron dynamics can be extended to include the spin degrees 
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of freedom, both in the hnear and nonhnear regimes. An outstanding ques- 
tion concerns the demagnetization processes observed in ferromagnetic thin 
films irradiated with femtosecond laser pulses, for which a clear theoretical 
explanation is still lacking. 

The field of optical control of spins in semiconductor nanostructures is also 
a very active research area. It is nowadays possible to fabricate and optically 
probe individual semiconductor quantum dots doped with one or more mag- 
netic impurities [94]. One of the major interest of this type of structure is the 
possibility to control magnetism via optical processes acting on the charge 
carriers. Thus, ferromagnetism becomes optically manipulable on an ultrafast 
timescale. This is particulary interesting for the elaboration of future fast- 
access magnetic storage devices. We are currently working on quasi one- and 
two-dimensional nonparabolic quantum dots containing up to four electrons 
and doped with a finite number of localized magnetic impurities. Within the 
framework of the CI method and the Anderson model, we aim at investigating 
the influence of the impurities on the energy spectra and oscillator strengths 
with special emphasis on the breakdown of the Kohn theorem. 

Finally, another procedure that has attracted particular attention over the 
last decade is the low-density doping of semiconductor nanostructures with 
magnetic impurities such as manganese ions. The resulting materials (named 
DMS, for diluted magnetic semiconductors) can display Curie temperatures 
as high as 80K [95], and possibly larger [96]. The spin of the Mn ions is coupled 
to the spin degrees of freedom of the electrons and holes, whose dynamics can 
be optically excited. DMS thus offer the possibility of using laser pulses to 
control the magnetization dynamics of semiconductor nanostructures. 

Given the wealth of fundamental issues and practical applications, the 
interplay of charge and spin effects in nanosized objects is bound to remain a 
major area of research in the coming years. 
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